Chirality distribution and transition energies of carbon nanotubes 
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From resonant Raman scattering on isolated nanotubes we obtained the optical transition energies, 
the radial breathing mode frequency and Raman intensity of both metallic and semiconducting 
tubes. We unambiguously assigned the chiral index (711,712) of ~ 50 nanotubes based solely on a 
third-neighbor tight-binding Kataura plot and find wrbm = (214. 4±2) cm -1 nm/d-|-(18.7±2) cm -1 . 
In contrast to luminescence experiments we observe all chiralities including zig-zag tubes. The 
Raman intensities have a systematic chiral-angle dependence confirming recent ab-initio calculations. 

PACS numbers: 78.67.Ch, 73.22.-f, 78.30.Na 



The successful preparation of single- walled carbon nan- 
otubes in solution where the tubes are prevented from 
rcbundling has opened a new direction in carbon nan- 
otube research Q, Q, S B ■ Strong luminescence by di- 
rect recombination from the band gap was detected in 
these isolated tubes, whereas in nanotubc bundles no lu- 
minescence is observed. The electronic structure of car- 
bon nanotubes and the optical transition energies vary 
strongly with their chiral index (711,712) 0|- Because the 
synthesis of nanotubes with a predefined chiral index has 
not been achieved so far, luminescence experiments were 
carried out on tube ensembles with unknown composition 
of chiral angles. Several attempts to assign the chiral in- 
dex (m, ti 2 ) to the experimentally observed luminescence 
peaks were reported 0, El ■ With a unique assign- 
ment, one could validate and possibly revise theoretical 
models of the electronic band structure. Moreover, such 
an assignment would allow to characterize the tubes after 
their production and to control their separation 

Bachilo et al. suggested an (711,712) assignment of the 
first and second transition energies in semiconducting 
tubes 0. Their assignment is based on pattern recogni- 
tion between experiment and theory in a plot of the sec- 
ond transition (excitation energy) versus the first transi- 
tion (emission energy) 0. The patterns, however, were 
not unique, and the frequency of the radial breathing 
mode (RBM) was used to find an anchoring clement that 
singles out one of the assignments. Surprisingly, zig-zag 
tubes were not detected in these luminescence experi- 
ments. Bachilo et al. concluded that the concentration 
of tubes with chiral angles close to the zig-zag direction 
was very low in the sample ■ 

The electronic transition energies of metallic nan- 
otubes cannot be detected by luminescence experiments. 
An elegant approach is to record Raman resonance pro- 
files 0,0, El, El, with maximum intensity close to the 
real transitions in the electronic band structure. Reso- 
nance profiles from nanotubes in solution were first re- 
ported by Strano et al. |l4j : their (711,772) assignment to 
the transition energies was based on the RBM frequency 



to tube diameter relationship of Ref . • The resonance 
profiles of the so-assigned RBM peaks were then used to 
find an empirical expression for the transition energies in 
metallic tubes. 

In this paper we present the transition energies of 
both metallic and semiconducting nanotubes by reso- 
nant Raman spectroscopy. Plotting the resonance max- 
ima as a function of inverse RBM frequency, we obtain 
an (711,712) assignment without any additional assump- 
tions. From our assignment we fit c\ = 214.4 cm" 1 !!!!! 
and C2 = 18.7cm -1 for the relation between diameter 
and RBM frequency. We observed several semiconduct- 
ing tubes that were not detectable by luminescence. Our 
results show that the electron-phonon coupling strength 
increases systematically for smaller chiral angles. Con- 
clusions about the distribution of chiral angles in a sam- 
ple based solely on luminescence intensity lead to incor- 
rect results; in particular, zig-zag tubes are present in 
nanotube ensembles. 

We performed Raman spectroscopy on HiPCO nan- 
otubes with diameters d w 0.7— 1.2 nm [J^. We dispersed 
the tubes in D2O containing a surfactant, see Ref. |j] and 
used an Ar-Kr laser between 2.18 and 2.62 eV and two 
tunable lasers (1.85-2. 15eV and 1.51-1. 75eV). The spec- 
tra were collected with a Dilor XY800 spectrometer in 
backscattering geometry at room temperature. To ob- 
tain the Raman cross section from the measured inten- 
sity we normalized the spectra to CaF2 and BaF2 mea- 
surements taken under the same experimental conditions 
(integration time, laser power). This also corrects for 
the spectrometer sensitivity and the u> 4 dependence of 
the Raman process. The Raman susceptibility was cal- 
culated from the normalized spectra by dividing by the 
Bose-Einstein occupation number and the inverse phonon 
frequency 0|. The latter was omitted in Fig.^(a) for a 
better representation. 

In Fig. H( a ) we show a contour plot of all Raman 
spectra, i.e., the Raman scattering power as a function 
of inverse RBM frequency (1/wrbm) and excitation en- 
ergy. When tuning the excitation energy, the RBM peaks 
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FIG. 1: (Color online) (a) Contour plot of the Raman cross section of the RBM as a function of excitation energy and reciprocal 
RBM frequency. In (a) and (b) the dotted and dashed lines connect maxima originating from tubes of the same branch. In 
each branch the member with the largest RBM frequency is labeled, (b) Kataura plot from experimental results (filled circles) 
and third-neighbor tight-binding calculations (open circles). Gray and black indicate semiconducting (2nd and 3rd transitions) 
and metallic tubes, respectively. The plot consists of theoretical (right, top axes) and experimental data (left, bottom axes). 
The error in experimental transition energies is mostly smaller than 0.02 eV. For the (9,1) branch, the vertical arrows indicate 
for each member the assignment to its theoretically predicted point. The increasing softening compared to the theoretical 
transition energies with smaller chiral angle is clearly seen, (c) Integrated Raman intensity as a function of excitation energy 
for the tubes in the (12,1) branch. 



appear and disappear in groups of close- by frequencies. 
These are indicated by the dotted lines, where the later- 
assigned chiral index of the largest RBM frequency of 
each group is given. To obtain the optical transition en- 
ergies, we use the resonance profiles, see, e.g., Fig. [2(c) 
for the group beginning with the (12,1) tube. The lines 
are fits to the first-order Raman cross section includ- 
ing incoming and outgoing resonance. From these fits 
we obtain directly the transition energies, which are by 
ps lOmeV (0.5wrbm) smaller than the energy of the res- 
onance maximum. For decreasing RBM frequency, the 
resonance energy first increases and then decreases again 
[dashed line in Fig. [2(c)] . The maximum intensity de- 
creases with decreasing RBM frequency. The described 
systcmatics are valid only for RBMs within a given group 
of tubes defined by the dotted lines in Fig. [2(a)- 

In Fig. [2(b) we plot the experimental transition en- 
ergies as a function of 1/o>rbm (solid circles, bottom 
and left axes). Since 1/wrbm is approximately propor- 
tional to the tube diameter, we thereby obtain an ex- 
perimental Kataura plot, which we compare to third- 
nearest neighbor tight-binding calculations Q (open cir- 
cles, upper and right axes). For semiconducting tubes, 
the measured transition energies are in excellent agree- 
ment with the E22 energies measured by Bachilo et 
al. Q. From a comparison between the experimental 
and the calculated Kataura plot we were able to assign 
the tube chiralities in the following way. Since we did 



not want to make any assumptions about the coefficients 
ci and C2 in the diameter-RBM frequency relationship, 
wrem = (ci/d+ C2), we varied c\ and C2 to find the 
best match between experimental and calculated data in 
Fig. ^(b). Varying c\ and C2, respectively, corresponds 
to stretching and displacing the theoretical Kataura plot 
along the 1/wrbm axis. In addition we displaced the plot 
on the energy scale. The origin of this overall energy off- 
set is, among others, found in the tight-binding param- 
eters, which were fitted to LDA calculations believed to 
underestimate the transition energies by typically 10 %. 

Figure [2(b) shows the best (72.1,72.2) assignment of 
experimental transition energies and RBM frequencies. 
The dashed lines in the theoretical Kataura plot indicate 
branches that are formed by tubes with 2ni + 72-2 con- 
stant. The indices of neighboring tubes in a given branch 
are related through (rii + 1, 712 — 2) or (tji — 1, n-i + 2), 
as long as 72. 1 > n^. For example, the branch starting 
with the (11,0) tube contains from small to large diam- 
eter (11,0), (10,2), (9,4), and (8,6), and belongs to the 
v = (ni — ri2)mod3 = —1 family. The chiral angle 6 of 
the tubes within a branch increases towards larger diam- 
eter. Zig-zag tubes (6 = 0°) arc always at the outermost 
position of a branch. The RBM groups in Fig. [dot- 
ted lines] directly correspond to these branches, see, e.g., 
the one beginning with the (9,1) tube. 

We find a very good match between the patterns of 
the calculated diameter and of the experimental values 



3 



TABLE I: Measured Wrbm and E22 for the branch of the 
(11,0) tube. See also supplementary material [tHI . 



chiral index 


(11,0) 


(10,2) 


(9,4) 


(8,6) 


wrbm (cm -1 ) 


266.7 


264.6 


257.5 


246.4 


E22 (eV) 


1.657 


1.690 


1.72 


1.73 



1/^RBM- The experimental transition energies, how- 
ever, deviate systematically from the calculations for the 
branches on the low-energy side of the semiconducting 
and of the metallic transitions. This deviation increases 
for smaller chiral angles. In other words, the experimen- 
tal transition energies bend downwards from the calcu- 
lated transitions with decreasing chiral angle. Such a 
strong deviation between theory and experiment is not 
seen for the upper Kataura branches, where the tubes 
belong to the v = +1 family. This was observed in lumi- 
nescence experiments as well and will be discussed below. 

We considered alternative assignments by displacing 
the experimental data along the sets of semiconducting 
and metallic tubes in Fig. H(b) to the left and to the 
right. None of them yields a good agreement regarding 
the data points within the (n\ — 1, n,2 + 2) branches: Ei- 
ther the number of tubes differs between the theoretical 
and experimental branch or some points are eminently 
displaced horizontally from the theoretical ones. Having 
found the zig-zag tubes we can exclude these alternative 
assignments. Table|I]summarizes the assignment of chiral 
indices to the measured RBM frequencies and transition 
energies for the branch of the (11,0) tube [l7| . 

From Fig.^(b) we can now fit the coefficients of the lin- 
ear relationship between wrbm and l/d of the assigned 
nanotubes. The linear fit based on 45 identified tubes 
yields wrbm = c-i/d + C2 with c\ = (214.4 ± 2) cm nm 
and C2 = (18.7 ± 2) cm -1 . Here, the tube diameter is 
geometrically determined by d = ao^n^ + ri\ri2 + nJ/71", 
using a graphite lattice constant ao = 2.461 A. c\ and C2 
differ somewhat from the coefficients found in Ref. Q for 
the same type of samples, because we used many more 
RBM frequencies for the fit and calculated the diameter 
from a smaller a n. T he coefficients are similar to theoret- 
ical predictio ns ll 3. Il9i| but different from other experi- 
mental work |l2l l20j. In Ref. |l2j . c\ and C2 were used 
as free parameters to find the assignment; moreover, the 
experimental information on the transition energies was 
not included. In Ref. j2lj, the chiral index assignment 
also depends on c\ (C2 = 0). Moreover, the tubes are 
on a substrate, which might alter the RBM frequencies. 
In contrast, our assignment is not based on a choice for 
Ci and C2- For the first time, they are obtained by a lin- 
ear fit after the assignment was performed. For example, 
the (11,0) tube has an RBM frequency of 266. 7cm^ 1 and 
E22 = 1.657 eV regardless of its exact diameter or any fit- 
ting procedure, see Table [I] 



In addition to the semiconducting tubes, we directly 
obtained the transition energies of metallic tubes. Our 
assignment of the metallic tubes to RBM frequencies 
agrees well with Strano et al. |l4l. The empirically ob- 
tained transition energies in Ref. 14], however, underes- 
timate the experimental values in Fig. ^(b). Moreover, 
our data show a stronger bending of the branches towards 
small chiral angles. This discrepancy comes mainly from 
the presence of pairs of close-by transition energies in 
chiral metallic tubes [2]]]. This ambiguity led to an in- 
correct assignment of single data points to the upper or 
lower transition. 

Given the unique (ni, 77,2) assignment of the RBM and 
the resonance energies, we can now examine the chirality 
dependence of the transition energies and of the RBM 
intensity. The optical transition energies of carbon nan- 
otubes are roughly proportional to l/d 0, 0, |U H3 • 
Tight-binding calculations predict deviations from a pure 
l/d dependence as a function of chiral angle I2J. They 
lead to the branches in the Kataura plot and are more 
pronounced if third-nearest neighbors are included [{J. 
Still, the band gap energies of tubes with small chiral 
angle (zig-zag tubes) are systematically lower in first- 
principles calculations than in zone folding |24j|. The 
branches in the experimental Kataura plot of Fig. H(b) 
are an experimental verification of the predictions from 
ab-initio calculations. 

The chiral-angle dependent softening of the transition 
energies with respect to the third-order tight-binding 
calculations is due to rchybridization of the it and a 
bands 0, 0, |2f| . It is stronger for the states originat- 
ing from between the K and M point of the graphite 
Brillouin zone in the zone-folding approach and weaker 
for states from the other side of the K point ji^. In 
semiconducting tubes, the value of v = [n\ — 71,2) mod 3 
determines from which side of the K point the electronic 
states originate for a given optical transition. The tubes 
in the lower branches of the E22 transitions in semicon- 
ducting tubes, i.e., of the branches beginning with (9,1), 
(11,0), and (12,1) in Fig. [2(b), have v = -1. The soft- 
ening for small-chiral-anglc tubes in these branches com- 
pared to the tight-binding value is quite strong. The 
tubes in the upper part of the same set of transitions 
have v — +1 and are less affected by the softening. For 
example, the experimental transition energy of the (10,0) 
tube [8 — 0°, v = +1] matches the theoretical value very 
well [Fig. n](b)]. Taking this softening of the transition 
energies in the —1 families into account, the agreement 
of our data with the theoretical predictions is excellent. 

In general, the RBM signal was strong for nanotubes 
with v = — 1 and also from the lower branches of the 
metallic transitions in Fig.^(b). The intensities were by 
a factor of four to ten weaker for tubes with v = +1. 
This observation confirms ab-initio calculations of the 
electron-phonon coupling that predicted the magnitude 
of the electron-phonon matrix element to alternate with 
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FIG. 2: (Color online) Raman intensity as a function of the 
chiral angle for three nanotube branches (solid circles) with 
(711,712) as indicated. Open circles: calculated Raman inten- 
sity, see text, (a) and (b) contain semiconducting, (c) contains 
metallic tubes. 

v j2|| . The relative Raman intensity of the RBM can 
also be used to discriminate between the two families of 
semiconducting tubes. 

The assignment of the semiconducting tubes in 
Fig. H](b) corresponds to the one found by Bachilo et 
al. from luminescence experiments 0- In contrast to 
the luminescence results, which reported a maximum in- 
tensity for closc-to-armchair tubes and no emission from 
zig-zag tubes, we clearly observed zig-zag or close-to-zig- 
zag tubes as well. The (13,0), (11,0), and the (10,0) 
tube show that zig-zag tubes are present in the sample. 
These tubes as well as the (14,1) tube were not observed 
by photoluminescence. An important conclusion from 
our work is that the absence of photoluminescence from 
these tubes thus does not imply the preferential growth 
of armchair tubes. 

In Fig. |21 we show the resonance maxima for tubes be- 
longing to the same {n\ — 1, ?i2 + 2) branch. The Raman 
intensity increases with decreasing chiral angle 9 and is at 
maximum for 9 « 10— 15°, except for (a). Theoretically, 
the Raman amplitude is proportional to the electron- 
phonon coupling times the square of the optical absorp- 
tion strength [lfij . Machon et al. found that the 
electron-phonon coupling of the RBM decreases strongly 
with increasing 9. We therefore model the electron- 
phonon interaction by a linear function of the chiral angle 
with a three times stronger coupling for zig-zag than for 
armchair tubes as found in Ref. |26| . We approximate 
the optical absorption strength of a tube by its experi- 
mental photoluminescence intensity 0, i.e., we assume 
the absorption and 

emission probability to be the same. The relative Ra- 
man intensities calculated with this model are in good 
agreement with experiment (open dots in Fig. [2J. The 
Raman response of zig-zag tubes is enhanced compared 
to their luminescence intensity because of their strong 
electron-phonon coupling. On the other hand, the Ra- 
man intensity of zig-zag tubes is smaller than for 9 « 10° 
due to the small absorption coefficient. Thus, our data 
are completely consistent with a uniform chirality distri- 



bution in the sample. 

In conclusion, we assigned the chiral indices to s=s 50 
measured RBM frequencies and transition energies by 
resonant Raman spectroscopy. In contrast to all pre- 
vious work our assignment is independent of the co- 
efficients ci and C2, which we fit only after assigning 
the chiral index to a particular RBM. The largest Ra- 
man intensity was measured for tubes with chiral angles 
around 15° or smaller, which is in agreement with the- 
oretical predictions and implies that the chiralities are 
evenly distributed. Moreover, our results confirm that 
the RBM intensity in semiconducting tubes depends on 
the (rii — 712) mod 3 family. The transition energies de- 
viate from zone-folding predictions with decreasing chi- 
ral angle, which, in particular for metallic tubes, was 
strongly underestimated in earlier work. 
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heimer Fund and Newnham College. 
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